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Abstract 



Recently the ground state and some excited states of the half-filled case of the 
Id Hubbard model were discussed for an open chain with L sites. Authors con- 
sidered the case when the boundary site has a negative chemical potential -p and 
the Hubbard coupling U is positive. They have shown by an analytic method 
that when p is larger than the transfer integral some of the ground-state solu- 
tions of the Bethe ansatz equations become complex-valued. They have found 
that there is a surface phase transition at some critical value p c ; when p < p c all 
the charge excitations have the gap for this case, while there exists a massless 
charge mode when p > p c . To find whether this "surface phase transition" 
is of the first order or of the second order we have used the entanglement en- 
tropy concept. The entropy and its derivative has a discontinuity there, so this 
transition is of the first order. 
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1 Introduction 



Quantum critical phenomena are phenomena where the transition is driven by 
quantum effects in the system of many electrons and bosons. Quantum phase 
transition with nonstandard behaviour was described in [1], it is a phase transi- 
tion and a critical behaviour deviating from the standard scheme based on the 
mean field theory renormalising only the mass of the critical excitations and 
with a perturbative scaling renormalisation of the interaction strength. 

Criticality in ld-Hubbard model was studied recently using entanglement 
scaling [2] . Entanglement is a quantity which characterises nonlocal correlations 
in quantum systems. There is a connection between the entanglement of a 
many particle system and the quantum critical phenomena. Discontinuity in 
the entanglement derivative seems to be [2] associated with the first order QPT, 
singularity with the second order. Entanglement describes how the associated 
nonlocal quantum correlations influence the critical behaviour of a quantum 
phase transition. As noted by the authors [2] to what extent their results can 
be generalised to other quantum systems is yet to be answered. In this paper 
we will discus entanglement for the ld-Hubbard model with the boundary site 
with a negative chemical potential -p and the Hubbard coupling U positive for 
a finite number of atoms in the chain. In general we have found that there is a 
1-st order "surface phase transition" phenomena present. 

2 Quantum Critical Phenomena and Entangle- 
ment. 

The Hubbard model displays many features which are close to real properties of 
materials with electron-electron interactions. Some rigorous results concerning 
this model are described in [3]. These results and the entanglement concept 
mentioned in the preceding section enable to study the quantum critical phe- 
nomena in ld-Hubbard model. Let us describe them shortly. 

Exact solution of the Hubbard model with boundary hoppings and fields 
was found in [4]. The Bethe ansatz method was used. It was found that for 
certain values of the boundary hoping integrals and fields the ground state 
contains boundary bound states. While in [2] the authors assume periodic 
boundary conditions, in [4] the boundary conditions are such that periodicity is 
broken. The Hamiltonian in [2] is transnational invariant and conserves particle 
number and the z component of the total spin. When the boundary fields and 
hoppings are present then this symmetry is broken. However boundary fields 
and hoppings are expected to play no role in the limit of large lattice, L — > oo. 
The reduced density matrix pA for a single site A is diagonal in the chosen 
basis for periodic boundary conditions. For boundary fields and hoppings or 
other boundary effects present we expect that it is diagonal in the limit of large 
lattice, L — > oo. 
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The single-site entanglement entropy of the ground state \ipo > can be writ- 
ten [2]: 



E = -w log 2 {w ) - w u log 2 (w u ) - w d log 2 {w d ) - w 2 log 2 {w 2 ) 



(1) 



where 



w 2 =< njuUjd > 



(2) 



W a =< Hja > — w 2 , 



WQ = 1 - W u - Wd - W 2 . 



where a = u,d are spin up and down indices. 

In the limit of large |u| = ||4| >> the Hellman-Feynman theorem can be 
used, < ^ >o= ^f-- The Bethe ansatz result for the ground state energy gives 
energy proportional to L, the number of sites in the lattice. There is present 
also a contribution to this energy from the surface which is of the order of L°, 
and contributions of the order of L _1 , see in [5]. 

3 Hubbard Model in Id with the Boundary Site 
with a Negative Chemical Potential -p and the 
Hubbard Coupling U Positive 

The surface energy for the boundary fields was calculated exactly for the XXZ 
chain by the Bethe ansatz in [6]. The energy is an analytic function of the field, 
where h = hi = —h 2 for the case in which there are opposite oriented fields h\ 
and —h 2 on boundaries. While there are identified "critical" fields h\ c and h 2c 
as fields below which and above which respectively (in absolute values) there are 
present 1-string boundary states in the chain, for fields with values in between 
hie and h 2c there are no boundary strings. Then using the entanglement we 
find that there is no discontinuity or divergence in its derivative for this case. 
So there is no quantum phase transition changing boundary fields and thus 
the ground states between the boundary-string state and the state without the 
boundary string state. The XXZ model is related to the Hubbard model by a 
canonical transformation. Thus we can expect that at the spin sector of the 
solutions of the Hubbard model with boundary fields and hoppings there is no 
" surface phase transition" . 

Let us now concentrate on the charge sector. In [7] the Hubbard model is 
discussed for the ground state and some excited states of the half-filled case for 
an open chain with L sites. Only one of the boundary sites has a different value 
of chemical potential in this case. Authors [7] considered the case when the 
boundary site has a negative chemical potential -p and the Hubbard coupling 
U is positive. They have shown by an analytic method that when p is larger 
than the transfer integral some of the ground-state solutions of the Bethe ansatz 
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equations become complex- valued. They have found, among other findings, that 
there is a surface phase transition at some critical value p3 C ; when p < p-i C all 
the charge excitations have the gap for this case, while there exists a massless 
charge mode when p > p^ c . The authors did not discussed whether the surface 
phase transition is of the first or of the second order in the sense discussed above 
for the entanglement. 

To find whether this "surface phase transition" is of the first order or of the 
second order we can use the entanglement entropy described above. We have 
found using the results of [7] that this entropy has a discontinuity there for 
zero magnetic field. Consequently the difference of derivatives of entanglement 
entropy with respect to the chemical potential is nonzero at p = ps c calculating 
it in the corresponding limiting procedure from below and from above p = p^ c . 
This difference of derivatives is proportional to a (charge) susceptibility and 
to a coefficient which is finite with exception of two points and | for w 2 . 
Thus in general there is a 1-st order "surface phase transition" phenomena 
present. The value of the quantity u>2 for p — P3c is in general different from the 
values corresponding to the mentioned two points for larger values of u when 
correlations of fluctuations of the number of electrons with spin up and with 
spin down is nonzero. For a finite chain with L finite we have for p > pz c a hole 
at the boundary site 1, and for p < p 3c there is an electron at this boundary site 
1 [7]. As we can see there is no phase transition in the thermodynamic limit. 
That is the reason why we speak about the "surface phase transition" in this 
model. Corrections to entanglement entropy are not logarithmic here, however 
again reflect a change in the number of local states accessible to the system at 
the transition as it was described above. 
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